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We explore the energy dynamics of a two dimensional (2D) harmonically trapped Bose-Einstein 
condensate inside a box potential, excited by a moving red-detuned laser potential (RDLP). For an 
RDLP velocity v less than a critical value Vc, energy oscillations are observed to begin simultaneously 
with the motion of the RDLP. For v > Vc, these oscillations are delayed through a transient in the 
energy. At the end of the delay time tons, the energy oscillations are regenerated again, and tons is 
found to depend on v through universal relations for two cases: one for an RDLP depth sufficient to 
break off a BEC fragment, another for a depth insufficient to split the BEC. In the case of splitting, 
tons exactly equals the time it takes the BEC fragment (dragged by the RDLP trough) to reach the 
hard wall of the box potential. 



a. Introduction: The investigation of a Bose-Einstein 
condensate (BEC) excited by an obstacle produced by 
a laser beam has received considerable theoretical and 
experimental interest in the last few years [Illl6j. So 
far, mostly repulsive obstacles have been created experi- 
mentally by a blue-detuned laser beam [HHS]- Recently, 
however, a new method has been suggested to create at- 
tractive obstacles using a red-detuned laser beam instead 
[17j . In the latter, vortex dipole generation was shown to 
be possible with two red-detuned Gaussian beams. The 
success of using a red-detuned laser for the nucleation of 
vortices as in Aioi et al. [17], and the lack of literature 
therein, provided the main motivation for this work. In 
another quest, we were motivated to search for a uni- 
versal behavior (UB) in our system, as UBs have been 
demonstrated earlier in trapped BECs [HI fT8H22] . For 
example, Raman et al. fl5' revealed a universal behavior 
connecting the thermal fraction of the condensate with 
the velocity of the laser. We then set out to search for 
a UB connecting the onset time for solitonic nucleation 
and the velocity of a moving attractive obstacle created 
by a red-detuned laser potential (RDLP). 

Our chief goal in this Letter, then, is to demonstrate 
the presence of two UBs connecting the "onset time" tons 
for solitonic excitations with the velocity of the moving 
laser beam, both for a shallow and an RDLP. The tons 
is recorded in the transient stage of the simulation [53], 
from the start of laser motion until the energy begins to 
display an oscillatory behavior. As we have outlined ear- 
lier 23^, and in correspondence to Parker et al. [21], the 
energy of solitons oscillates with time. In addition, we 
show that a deep RDLP suppresses soliton excitations. 
To the best of our knowledge, our current findings per- 
taining to the energy dynamics in a trapped BEC excited 
by an RDLP have not been reported earlier. 

In our previous publication [23] . we have demon- 
strated that a blue-detuned laser potential (BDLP) 
moving inside a harmonic trap plus hard-wall box- 
potential (HWBP) boundaries can lead to fascinating 



self-interfering matter wave patterns. In this work, we 
invert the "sign" of the laser potential (LP) to make it 
attractive and investigate the dynamics along the same 
lines as in our Ref . ^23j . We follow the method of splitting 
and merging a number of BEC fragments [T4 1I25H28] . For 
example, Scherer et al. [TJ] successfully observed vortex 
formation by splitting a harmonically trapped BEC into 
three sections using a blue-detuned laser, and then al- 
lowing those sections to merge. Chang et al. [28] has 
demonstrated experimentally, using a BDLP, that soli- 
tons and dispersive shock waves (DSWs) can be produced 
by merging and splitting of BECs, respectively. DWs 
were also produced by merging BECs using a RDLP. By 
using a "deep" RDLP well, we separate off a BEC frag- 
ment from the main BEC and move it toward the HWBP 
boundary. Once it is reflected from the hard wall, after 
the RDLP leaves the trap, it merges again with the cen- 
tral part of the BEC. In essence, another goal of this 
Letter is to motivate experiments which explore the dy- 
namics of a trapped BEC excited by an RDLP. 

Our key results are (i) tons depends on the RDLP 
velocity v only, and is independent of RDLP width and 
depth and independent of the coupling strength M be- 
tween the bosons, i.e., it displays universal behavior. 
Two universal behaviors for the tons vs v were found as 
follows: a) for an RDLP well deep enough to split the 
BEC, tons was found to be exactly equal to the time it 
takes the separated fragment to reach the HWBP bound- 
ary from the BEC center, that is tons = 10/w; b) for 
a shallow RDLP, insufhcient to split the BEC, a dif- 
ferent universal behavior was obtained for tons, namely 
tons — 3.3/?;. We verified [i(a)] by fitting a classical 
curve for the relationship between tons and v, which is 
tons = d/v, where d = 10 is the distance from the center 
of the trap to the HWBP boundary, (ii) An RDLP sup- 
presses solitonic excitations, since before the separated 
BEC fragment reaches the HWBP boundary, no oscilla- 
tions appear in the energy dynamics. 
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FIG. 1: Two dimensional (2D) density \ip(x,y;t)\'^ at t = 4.5. 
The system is a hard core Bose gas in a 2D harmonic trap 
plus box-potential boundaries [23] excited by a red-detuned 
laser potential with A/" = 10, A = —30, P — 4, v — 2, and 
K — 1. A is in units of htOho, P in {aho)~^ , v in Uho, t = TUJho 
is unitless, Af is in units of {\/2af^^)~^ , and \ip{x,y;t)\^ is in 
{aho)~'^- X and y are in units of aho- 



b. Method: The model we use is the same as in our 
previous work [23 , except that now the LP is a trough in- 
stead of a barrier. The split-step Crank Nicolson method 

is applied to solve the 2D time-dependent Gross- 
Pitaevskii equation (TDGPE) [Eq.{l) in Ref.[23]]. The 
energy E{t) (function of time) is computed using Eq.(7) 
in Ref.[23!. The coupling constant TV, LP depth A, LP 
velocity v, inverse width of the LP /3, total energy of the 
system and time tons, all have the same units as be- 
fore P^^: The lengths and energies are in units of the 
trap Qho — \/^/ {'2.mujho) and Tiujho, respectively. Hence, 
A is in units of fiuJuo, P in (o/io)"^, v in auo, t — nuho 
is unitless, JV is in units of (v^a^^)"^, and the density 
\ip{x,y;t)\^ is in (a^o)^^. We also use a spherically sym- 
metric harmonic trap. The atoms are attracted to the 
well of the RDLP, thereby creating a density peak as the 
atoms fill the trough. As the RDLP well moves, and for 
large enough \A\, larger than a critical value Ac > 0, it 
separates off a fragment of the BEG and carries it away 
from the main part. As before, our system is cut off by 
HWBP boundaries [23], but the RDLP can continue its 
motion and leave the boundary. As the fragment reaches 
this boundary, it bounces off and returns in the opposite 
direction. 

c. Results and Discussion: Figures [T] and |2] display 
the 2D BEG densities at different times t = 4.5 and 
5.4, respectively. The system is a trapped BEG in an 
isotropic 2D harmonic trap with HWBP boundaries, as 
in Ref. [23], with interaction strength TV — 10, potential 
depth A ~ —30, inverse width f3 — A, velocity v — 2, 
and K — 1. In Fig. [TJ the separated BEG fragment 
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FIG. 2: As in Fig. [TJ but for t = 5.4. The units are as in Fie. 

El 

dragged by the RDLP collides with the hard wall at 
y = 10 and t = 4.5. In Fig. [2] the BEG fragment has 
bounced off the wall and is on its way back to the center 
of the trap. The RDLP has aheady left the HWBP. Since 
the BEG fragment is no longer trapped by the RDLP, it 
increases in width and declines in amplitude as it merges 
into the central BEG. By monitoring the energy dynam- 
ics of this system in the transient stage of the simulation, 
one observes an oscillatory pattern indicative of soliton 
excitations. This is demonstrated in Fig. |3] which dis- 
plays the energy as a function of time, E{t), for a system 
with TV = 10, A = —30: (open triangles) v = 0.3, 
/3 — 2; (open squares) v — 0.4, (5 — 2; and (open cir- 
cles) V — 2, (3 — 4. These oscillations are similar to 
those in Figs. 11-13 of our previous article 23 , where a 
BDLP was used. However, whereas in the case of a BDLP 
the energy begins to oscillate once the BDLP is moved, 
in the case of an RDLP the energy oscillations suffer a 
delay for a certain time if the velocity of the RDLP ex- 
ceeds a critical value Vc- Thus if we begin with a small 
velocity v — 0.3 and (3 — 2, the energy oscillations begin 
simultaneously with the motion of the RDLP, and a reg- 
ular oscillatory pattern is observed (open triangles). On 
increasing v slightly to 0.4 for the same (3, these oscilla- 
tions are damped significantly, and delayed from i = to 
t = tons = 22.5, after which they are reenergized (open 
squares). (The arrows in the figures indicate where the 
oscillations begin after they are delayed.) Gonsequently, 
then, there is a critical velocity, Vc, between v = 0.3 and 
V = 0.4 at which the energy oscillations begin to suffer 
a time delay. This Vc depends on the parameters of the 
system. However, the precise determination of would 
require a substantial number of runs, and we therefore 
suffice with an estimate of Vc- It is also possible, that 
the system for v = 0.4 is in a transitional state. In com- 
parison, for V — 2 and a different (3 — 4 (open circles), 
the oscillations begin at a time tons — 4.3 and decay 
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FIG. 3: Energy dynamics E{t) for the same systems as that 
in Fig. IT] but for v — 0.3, P = 2 (open triangles); v = 0.4, 
(3 = 2 (open squares); and w = 2, /3 = 4 (open circles). The 
arrow indicates the onset time for solitonic nucleation, tons- 
E{t) is in units of huJho and the rest of the units are as in Fig. 
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FIG. 4: As in Fig.js] but for v = 0.2, /3 = 4 (sohd triangles); 
V = 0.3, /3 = 4 (open squares); and v — 0.4, 13 — 4 (open 
circles). E{t) is in units of huho and the rest of the units are 
as in Fig.[l] 



afterwards for t ^ 10. The decay in these oscillations 
signals that upon a merger of the two separated BECs, 
the solitons become unstable and decay into vortex rings 
in the high density buldge formed during the merger |28j . 
Likewise, Raman et at |15| found that above a critical 
velocity the superflow around a BDLP becomes unstable 
against the creation of quantized vortex lines. They also 
reported evidence for the existence of a critical velocity, 
where the EEC gas enters into a new dissipative regime. 

For V > Vc, the initial energy of the system is neg- 
ative, because for the parameters chosen the negative 
potential-energy of the RDLP outweighs the combina- 
tion of kinetic, interaction, and trapping energies. As 



time progresses, E{t) increases and becomes positive, as 
the moving obstacle transfers energy to the system in 
the initial stages of the simulation [30]. Eventually, E{t) 
displays a transient similarly to the one reported earlier 
by Law et al. |30| for a repulsive potential. However, 
whereas in Ref. [30] the system is a homogeneous EEC, 
where the energy stabilizes after the transient, in our 
case the transient is followed by oscillations for the red- 
laser case. Thus, the principle feature in Fig. [3] above, 
is that a transient is also observed using an RDLP. In 
addition, according to Jackson et al. j5 , the creation of 
phonons or vortices increases the energy of the conden- 
sate, which plays a role in the rise of the energy in the 
transient regime. For a repulsive-obstacle velocity lower 
than a critical one Vc, energy is transferred to the system 
via phonon excitations; for a velocity larger than Vc, en- 
ergy is transferred via vortex formation. The oscillations 
in the energy following the end of the transient are, as 
has been outlined earlier [53] due to soliton/ vortex ex- 
citations. Therefore, an important finding here is that 
the RDLP suppresses the solitons until the RDLP leaves 
the HWEP, as manifested by the transient in the en- 
ergy. This is in stark contrast to the EDLP case which 
never suppresses solitons. After t > tons, the reflection 
of the separated EEC fragment off the HWEP boundary 
creates a shock wave in the system, which excites soli- 
tons, as the EEC fragment temporarily acts like a EDLP 
obstacle by itself, until it merges again with the cen- 
tral EEC part. Prior to this reflection, the RDLP pecu- 
liarly does not generate these excitations. The two EEC 
fragments "communicate" together by a shock wave, and 
cause these oscillations to arise. The backfiow around 
this EEC fragment is also a way in which the two BECs 
can communicate together [3T]. The effects of drag due 
to the interaction of the obstacle with the surrounding 
fluid have been explored earlier [5] ES] , and we shall 
explore these effects in a future extended version of this 
Letter. Next to this, when the RDLP leaves the system, 
the effective interaction jumps from being net attractive 
to net repulsive. This is, as if one "quenches" the Hamil- 
tonian "upwards" by a sudden reduction of the attractive 
forces. Due to the latter quench, the system produces 
topological defects. 

Figure |4] is the same as Fig. [3] except for a system 
with V — 0.2, /3 = 4 (open triangles), v = 0.3, (3 — A 
(open squares), and v — 0.4, /3 = 4 (open circles). Es- 
sentially, the same features are observed as in Fig. [3] a 
critical behavior between v = 0.2 and 0.3, and a reduc- 
tion in the "onset time" with an increase in v to 0.4. 

Next, we explore the possibility for the presence of 
universal behavior in the EEC, such as the ones reported 
earlier [TS ] IT8H22] . First, by changing v and keeping A 
and (3 fixed, we determined tons as a function of v (> v^) 
for A = —30, /3 = 4, and three different values of N dis- 
played in the top frame of Fig. [5] (solid squares) A/" = 0; 
(open circles) N = 10; (solid triangles) N = 20; (solid 
line) "fit" using classical formula f{v) = c?/w, where 
d = 10 is the distance between the center of the trap 
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once the BEC fragment dragged by the RDLP hits the 



FIG. 5: Onset time for solitonic nucleation tons versus tlie 
laser velocity v for the same system as in Fig. [l] with various 
parameters. Here A is deep enough to split the BEC into 
two fragments. Top frame is for fixed = — 30, /3 = 4 and 
difi^erent A/": (solid squares) — 0; (open circles) 10; (solid 
circles) 20; (solid line) fit f{v) = d/v. Middle frame is for 
fixed J\f = 10, /3 = 4 and different depths A: (solid squares) 
A = —20; (open circles) —30; (solid circles) —40. Bottom 
frame is for A = —30, A/" = 10 and different /3: (solid squares) 
P = 2; (open circles) 4; (solid circles) 6. Units are as in Fig. 
□ 



and the hard wall (see Figs. [T] and [2]). Second, the mid- 
dle frame of Fig. [5] is as in the top frame; but for differ- 
ent amplitudes A and fixed J\f — 10 and /? = 4: (solid 
squares) A = —20; (open circles) A = —30; (solid tri- 
angles) A — —40; (solid line) f{v). The bottom frame 
of Fig. [5] is for fixed Af = 10 and A = -30, but dif- 
ferent widths /3: (solid squares) /3 ~ 2; (open circles) 
/3 = 4; (sohd triangles) /3 = 6; (sohd hne) f{v). The 
basic feature in all these figures, is that for values of |^| 
high enough to separate off a BEC fragment as in Fig. 
[T the "onset time" for energy oscillations (solitonic nu- 
cleation) is independent of the parameters of the system. 
This demonstrates a universal behavior for the onset time 
as a function of the velocity v of the RDLP, which is al- 
most identical to the classical f{v) — d/v. Of course, 
it may seem obvious that the time it takes the RDLP to 
reach the HWBP boundary is simply d/v; but the sig- 
nificant feature here is that the energy oscillations begin 




FIG. 6: As in Fig. [5] but for various parameters A, j3, and 
as indicated in the figure. Here, A is not deep enough to split 
the BEC into two fragments. Units are as in Fig. |5] 



HWBP wall. This is also verified by the classical curve 
f{v). The peculiar feature in Fig. [s] bottom frame, is 
that even though for a lower /3 the width of the obstacle 
is larger, it takes the same time to initiate energy oscilla- 
tions as for a narrower one with larger /3. The peculiarity 
arises because it takes more time for a wide obstacle to 
"completely" leave the HWBP trap than a narrow obsta- 
cle, as measured relative to the HWBP wall. 

Finally, we display in Fig. [6] the behavior of tons with 
V for low enough \A\ < Ac, such that the RDLP is not 
strong enough to separate off a BEC fragment. The BEC 
dynamics in this case is totally different. Here \A\ was not 
deep enough to suppress solitons for times comparable to 
Fig. [5] By changing \A\, Af, and /3, we again observe (as 
labelled from cross to star) another universal behavior 
for tons, obeying tons — 3.3/w. Consequently, there are 
two universal behaviors for tons, one for \A\ > Ac at 
which a BEC fragment is separated, and one for < Ac 
for which the RDLP is not strong enough to separate a 
fragment! 

d. Conclusions: In summary, then, we have explored 
the energy dynamics of a harmonically trapped 2D BEC 
inside a HWBP excited by an RDLP. In general, it was 
found that an attractive obstacle leads to peculiar fea- 
tures in the energy dynamics, not observed earlier in the 
case of a repulsive obstacle [13] ■ It '^''^^ found, that the 
behavior of tons with the velocity of the RDLP v > Vc 
displays two universal behaviors according to whether 
the RDLP is deep or shallow. We believe that this re- 
search will provide a new direction in the investigation of 
trapped BECs excited by a red-laser potential, especially 
since there seems to be a lack in the theoretical and ex- 
perimental investigations applying the RDLP in trapped 
BECs. 
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